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Abstract 

Incompressible MHD turbulence is investigated under the presence of a uniform magnetic field Bo- Such a situation 
is described in the correlation space by a divergence relation which expresses the statistical conservation of the 
Elsasser energy flux through the inertial range. The ansatz is made that the development of anisotropy, observed 
when Bo is strong enough, implies a foliation of space correlation. A direct consequence is the possibility to derive 
a vectorial law for third-order Elsasser moments which is parametrized by the intensity of anisotropy. We use the 
so-called critical balance assumption to fix this parameter and find a unique expression. To cite this article: S. 
Galtier, C. R. Physique 11 (2010). 

Resume 

Moments d'Elsasser du troisieme ordre en turbulence MHD axisymetrique. La turbulence MHD in- 
compressible est etudiee en presence d'un champ magnetique uniforme Bo- Une telle situation est decrite dans 
l'espace des correlations par une relation de divergence qui exprime la conservation statistique du flux d'energie 
d'Elsasser a travers la zone inertielle. Nous faisons l'ansatz que l'anisotropie, observee quand Bo est suffisamment 
fort, implique un feuilletage de l'espace des correlations. Une consequence directe est la possibility d'obtenir une 
nouvelle loi vectorielle pour les moments d'Elsasser d'ordre trois qui est parametrisee par l'intensite de l'anisotro- 
pie. Nous utilisons l'hypothese d'equilibre critique pour fixer ce parametre et trouver une expression unique. Pour 
citer cet article : S. Galtier, C. R. Physique 11 (2010). 
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1. Introduction 

Despite its large number of applications such as climate, atmospherical flows or space plasmas, turbulence is 
still today one of the least understood phenomena in classical physics; for that reason any exact results appear 
extremely important [1]. The Kolmogorov's four-fifths (K41) law [2] is often considered as the most important 
result in three-dimensional (3D) homogeneous isotropic turbulence: it is an exact and nontrivial relation derived 
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from Navier-Stokes equations which implies the third-order longitudinal structure function. When isotropy is not 
assumed the primitive form of the K41 law is the divergence equation [3] 

-iv r -F HD (r)= £ , (1) 

where e is the mean energy dissipation rate per unit mass, r is the separation vector, F HD (r) = (SvSv 2 ) is 
associated to the energy flux vector and 5v = v(x + r) — v(x). Then, the K41 law may be seen as a non trivial 
consequence of equation (1) when isotropy is assumed; it is written as [2] 

-\er = {8vl), (2) 

where L means the longitudinal direction along r. Few extensions of such a result to other fluids have been made; 
it concerns e.g. scalar passively advected such as the temperature or a pollutant in the atmosphere [4] or space 
magnetized plasmas described in the framework of magnetohydrodynamics (MHD) [5], electron [6] and Hall [7] 
MHD. 

In this paper we investigate 3D homogeneous incompressible MHD turbulence for which the following divergence 
relation holds [5] 

-i Vr -F ± (r)= e ± , (3) 

where F±(r) = (<5z T (Sz+f) , z ± = v ± b are the Elsasser fields and e are the mean Elsasser energy dissipation 
rates per unit mass. When isotropy is assumed we obtain the exact law for 3D MHD [5] 

-*e ± r = {5zZ(8z ± ) 2 ), (4) 

which may reduce to expression (2) when the magnetic field is taken equal to zero. It is straightforward to 
demonstrate the compatibility between relations (3) and (4) by performing an integration of the former over a 
full sphere (ball). The same remark holds for the compatibility between expression (1) and the K41 law. 

To date the universal isotropic scaling relations discussed above have never been generalized to 3D homogeneous 
- non isotropic - turbulence (see however [8,9] for the latest progress). It is basically the goal of this paper to 
show that an exact relation may be derived in terms of Elsasser fields for axisymmetric MHD turbulence. This 
derivation is based on the ansatz that the space correlation is foliated when the field fluctuations are dominated 
by a uniform magnetic field. Note that a first analysis was made for such a problem in [8]. The main goal was 
the development of a tensorial analysis only for vectors v and b since the Elsasser fields, a mixture of a vector 
and a pseudo-vector, renders the study much more difficult. Then, the idea of foliation of space correlation was 
eventually introduced to derive a law for third-order correlations in v and b. In the present paper we show that 
the extension of the latter idea to Elsasser variables is possible - independently of their tensorial nature since we 
do not perform a tensorial analysis - and we derive the corresponding exact law. 



2. Impact of a mean magnetic field 

The influence of a large-scale magnetic field B on the nonlinear MHD dynamics has been widely discussed 
during the last fifteen years. The first heuristic picture of MHD turbulence proposed by Iroshnikov-Kraichnan 
[10,11] has been criticized and, nowadays, we know that under the presence of B we find turbulent fluctuations 
with larger fluctuating components in the direction transverse to Bo than along it, as well as different type of 
correlations along Bo and transverse to it [12,13,14,15,16,17,18]. In other words, the nonlinear transfer occurs 
differently according to the direction considered with a weaker non linear transfer along B than transverse to 
it, with possibly different power law energy spectra. An important concept introduced in the last years is the 
possible existence of a critical balance between the nonlinear eddy-turnover time and the Alfven time [19]. The 
former time may be associated to the distortion of wave packets whereas the latter may be seen as the duration 
of interaction between two counter-propagating Alfven wave packets. A direct consequence of the critical balance 
is the existence of a relationship (in the inertial range) between length-scales along (||) and transverse (_L) to the 
mean magnetic field direction (see also [20]). This relation, generally written in Fourier space, is 

fc„ ~ k 2 / 3 . (5) 
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In practice, numerical evidences of relation (5) may be found by looking at the parallel and perpendicular (to the 
mean magnetic field direction) intercepts of the surfaces of constant energy, either in physical space with second- 
order correlation functions [13,21] or in Fourier space with spectra [18]. Note that one generally takes a local 
definition for fcii by using the local mean magnetic field but it has been shown that a global definition (with the 
parallel direction along Bo) works quite well if Bo is strong enough [18]. Despite the limitation of direct numerical 
simulations a scaling relation between parallel and perpendicular length scales seems to emerge whose power law 
relation is compatible with the critical balance relation (5). Therefore, the idea of a general relationship between 
length scales during the nonlinear transfer (of energy) from large to small scales may be seen as a natural constrain 
for theoretical models. Basically, we translate this constrain as an ansatz for axisymmetric MHD turbulence which 
allows us to derive from equation (3) the equivalent of the four-fifths law. 

At this level of discussion, it is interesting to remark that the assumption of isotropy made to derive the exact 
law (4) is questionable in the sense that we never observe exactly isotropy. For example in [22,23] it was shown 
numerically that despite the absence of a uniform magnetic field (Bo = 0) deviations from isotropy are observed 
locally with the possibility to get a scaling relation between length-scales along and transverse to the local magnetic 
field. This local anisotropy is expected to be stronger at larger (magnetic) Reynolds numbers for which the exact 
law (4) is derived. Therefore, this exact law (4) should be seen as a first order description of MHD turbulence 
when Bo = 0. More precisely in the derivation of this law one should consider the decomposition 

F ± (r) = F± (r) + ,5F± ni (r), (6) 

where the first term in the RHS is the isotropic contribution to the vector third-order moment whereas the 
second term measures the deviation from isotropy. When the second term is of second order in importance then 
d^ini ^ ^iio an d the integration of relation (3) over a full sphere - with the application of the divergence theorem 
- gives the universal law (4). 

The derivation of a universal law from equation (3) in the general case of non isotropic turbulence is far from 
obvious. For example, one needs to find a volume V such that at its surface S the normal component F n of F is 
conserved. Then, one can perform an integration of equation (3) over this volume, apply the divergence theorem 
and obtain a simple expression independent of any parameter. In practice, that means one starts with 



UII^- F± ^ dv = s± III dv ' 



(7) 



which gives by the divergence theorem and after integration over the volume 

-\j J F ± (r)-dS = e ± V, (8) 
s 

and after projection on the surface vector dS 



-III 



Ft{v)dS^e ± V. (9) 



s 

If one assumes that F„(r) is constant on S then one obtains 



!*<•>// 



dS = -\F±(r)S = e ± V, (10) 



s 

which leads to the exact law 

Ft(r) = -4 £ ± ^. (11) 

The form (and even the existence) of such a volume V is still an open question. However, it is important to 
note that there exists an infinity of mathematical solutions of equation (3) but they depend on parameters which 
render the solutions non universal. For example we may have [24] 

F^r) = -4e(A ± pe p + (1 - 2A ± )ze*) , (12) 

where p and z are the cylindrical coordinates, and e p and e z are the corresponding unit vectors (with e z = Bo/Bo)- 
Note that the choice A ± = 1/2 gives the universal law for two-dimensional isotropic MHD turbulence, whereas 



3 



A = 1/3 leads to a radial vector and corresponds to the three-dimensional isotropic law [5]. Then, we may expect 
that relation (12) describes correctly anisotropic MHD turbulence when A G [1/3; 1/2] with stronger anisotropy 
when A ± is closer to 1/2. However, relation (12) does not satisfy the critical balance relation (5) for any values 
of A : indeed, for isotropic turbulence the energy flux vector is radial which may express the fact that energy 
cascades radially, whereas when a mean magnetic field is present it is not the case anymore and iso-contours of 
spectral energy are elongated in the perpendicular direction according to the power law (5) with an elongation 
more pronounced at small length scales (which means, in the correlation space, an elongation along the mean 
magnetic field direction). According to relation (12), we see that for a given distance r the energy flux ratio 
between a point along e z and another point along e p is equal to the following constant 

F ± (re z ) _ 1-2A ± 

F±(re p ) - A± ' {i6) 

This constant can be very small (when A ± is close to 1/2) but its precise value does not change the nature of the 
relation between these two fluxes which is linear. Therefore, it can only lead to a linear law dependence between the 
parallel and perpendicular intercepts of the surfaces of constant energy (the form of these surfaces being directly 
related to the intensity and direction of the energy flux). Note that if one considers a slightly different situation 
with points close to the e p and e z directions with energy fluxes F ± (re p + ee z ) and F ± (ee p + re z ) respectively 
(where e is a small parameter), the conclusion does not change drastically as long as r > e; when r becomes of 
the order of e then both energy flux vectors deviate significantly from the e p and e z directions which does not 
help for increasing anisotropy at small length scales which needs to have energy flux vectors preferentially along 
Bp. Expression (12) is the simplest solution among an infinity of axisymmetric solutions obtained by [24]. The 
expression that we shall derive here for the energy flux vector is another particular solution of this family which 
satisfies this time the critical balance assumption. 

In order to recover an anisotropic law of the type of (5) - which is a power law - it is necessary to reinforce 
the energy flux in the e p direction at small length scales. Then, the following statement is made that the energy 
flux vector has an orientation closer to the e p direction when the length scale decreases. This variation must have 
a power law dependence (with power law index n) in the length scale in order to be compatible with relation 
(5) which is also a power law. The value of n compatible with the index 2/3 in relation (5) may be determined 
with critical balance arguments (see Section 4). We will see that if we incorporate such a requirement in the 
analysis then we may derive a universal law in the sense that it does not depend on any (non physical) parameter. 
In practice, the energy flux vectors will belong to an axisymmetric surface S n in the three-dimensional space 
correlation (which means that F ± (r) is tangent to S„ for any points M' G S„; see Section 3 and Fig. 1). The 
manifold S„ is defined in such a way that the energy flux vectors tend to be perpendicular to e z when the distance 
separation goes to zero which means that turbulence tends to be bi-dimensional at small scales. As we will see in 
Section 6, the expected constant —2 for two-dimensional MHD turbulence is indeed recovered from the exact law 
when the small scale limit is taken. 



3. Foliation of space correlation 

From several theoretical and numerical analyses we know that MHD turbulence under the influence of Bo 
develops anisotropy that increases as the length scale decreases. Additionally, the rms fluctuations at a given 
separation distance r are more intense when r is perpendicular to B than when r is parallel to it. This property can 
be understood as a consequence of the critical balance relation (5) which provides a relationship between the length 
scales of the fluctuations parallel and perpendicular to the mean magnetic field. Following these considerations and 
those exposed at the end of Section 2, we make the ansatz that the energy flux vectors belong to two-dimensional 
surfaces S„ in the three-dimensional space correlation (which means that F (r) is tangent to S n for any points 
M' G S n ; see Fig. 1). Since the problem is axisymmetric, the manifolds S n must be of revolution about the (Mz) 
axis (with e z = B /-Bo; see Fig. 1). It is defined in such a way that the direction of F ± (r) tends to become 
perpendicular to e z when the distance separation r goes to zero. This variation of direction for F (r) should 
have a power law dependence in the length scale. Then, the axisymmetric manifold S n is defined by the following 
function 

z = MP) = Po(j^ • (14) 
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It is the simplest algebraic function satisfying the conditions fn(p) — > when p — > with a simple power law 
dependence between p and z. Without loss of generality we may already note that n must be greater than one to 
satisfy the anisotropic property (the energy flux vector getting perpendicular to B at small separation distance 
r). Finally, note that po is the value of p for which the angle between r and e z is n/4; therefore p/po may be seen 
as a way to delimit the correlation space into two domains where the direction of the separation vector r is closer 
to the transverse plane (xMy) or to the parallel direction e z (see Fig. 1). It is important to emphasize that the 




Figure 1. Left: We perform an integration of relation (3) over the manifold S n defined in the half upper space by the function 
fn(p) — Po(p/ Po) n with n > 1; note the use of the polar coordinates with r — (p, z). S n is a surface of revolution about the (M z) 
axis: on this Figure it appears as a "bowl" of axis of symmetry (Mz). The vector eT at point M' is tangent to the surface S n and 
perpendicular to the circle L n of radius p which has also (M z) for axis of symmetry. The curved vectors represent schematically 
the orientation of the energy flux (r) which flows towards the point M . For comparison the three-dimensional isotropic case is 
also represented (right) for which the energy flux flows radially. 

critical balance measured in MHD turbulence (with Bo > 0) is a situation towards which the nonlinear dynamics 
converges: it is the main state of the dynamics. In other words, deviations from this state may be found but are of 
second order in magnitude. In the same way, the assumption of a foliation of the space correlation (with relation 
(14)) means that one should write 

F ± (r) = F± 1 (r) + 5F± onfol (r), (15) 

where the first term in the RHS is the vector third-order moment which belongs to the foliated space correlation 
(see the schematic vectors in Fig. 1) whereas the second term corresponds to other vector contributions which are 
assumed (ansatz) of second order in importance, namely <5F* onfol <JC Ff ol . 

Equation (3) is integrated over the manifold S n of axis of symmetry (Mz). An illustration is given in Fig. 1 
where S n appears as a " bowl" . It gives 

-4e ± j JdS n = J J V r -¥ ± (r)dS n . (16) 

S n S„ 

By the Green's flux theorem (see Appendix) and after integration over the surface, we obtain 

-4e ± S n = j ¥ ± (r)-dL n , (17) 

circle 

where the line integral is performed along a circle L n of radius p and of axis of symmetry (Mz). On the example 
given in Fig. 1, it corresponds to the upper boundary of the "bowl". Note that dL n is an elementary vector which 
is normal to the circle L n and tangent to the surface S n (see Appendix). Then, one gets after projection 



5 



-4s ± S n = j F^(r)dL n , (18) 

circle 

where T means the tangent direction at point M' (see Fig. 1). The problem being axisymmetric, Fr(r) is unchanged 
along the circle L n of axis of symmetry (Mz); then we have 

-4s ± S n = F±(r) j dL n = F±(r)2np, (19) 

circle 

and thus 

-^=^W. (20) 
If we introduce the unit vector ei along the T-direction we obtain the vectorial relation 

- 2£ ^e T = F±(r), (21) 

TXp 

with 



_ e p + fn(p) e z _ e p + n (p/p»)" le z _ e p + wtange z ; , 
^1 + f k (pf ~ \/l + n2(p/p )2(«-i) ~ Vl + n2tan 2 e ' 



where 6 is the angle between r and the (xMy) plane (see Fig. 1). Note that for the foliated space correlation 
defined with relation (14) the general form of the divergence operator is 

v . ps l^) + l^ i (23) 
p al p o<p 

where <f> is the angle defined in cylindrical coordinates (note that by symmetry F$ — 0) and dT is the unit length 
along the tangent direction (see Fig. 1). The surface S n for a given p is defined as 



TTPo 



if""* 



S n = J 2npdT = J 2tv P ^1 + f{ l (pfdp = J 2npJl + ri^- 



X 

ny J y/l + X^dX, (24) 



n 2/(„ 





with 



2 / 2 /(n-l) 

The combination of the different expressions gives eventually the following vectorial law for Elsasser fields 



X = n 2 ^(j-)' = (™)" " = (ntan^) 2 ^"- 1 ) . (25) 



-2ffl £ ± pe T = F±(r), (26) 

where 

x 



7(X) = / Vl + X"- 1 dX. (27) 



4. Critical balance condition 



The vectorial relation (26) implies a parameter n that has to be determined. We shall fix n by a dimensional 
analysis based on the critical balance condition [19]. To investigate this idea we will restrict our analysis to the 
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inviscid, stationary MHD equations since basically we want an interpretation of relation (26) valid in the inertial 
range; we thus obtain 

z T • \7z ± = -VP* ±B a||Z ± , (28) 

where P* is the total pressure. By first noting that the divergence operator applied to (28) allows us to link the 
total pressure to the left hand side term, and second that z + ~ z~ for small cross-correlation; we then arrive to 
the nontrivial critical balance 

zfV r ~ Body , (29) 

which may also be written as 

where 6 is also the angle between the separation vector r and the (xMy) plane (see Fig. 1). As we see, relation 
(30) offers a direct evaluation of the r-direction: therefore, although the external magnetic field does not enter 
explicitly in the vectorial relation (26), it constrains - as expected - the direction along which the scaling law 
applies. If we now come back to relation (26), we may write (at first order for small length scales) the dimensional 
relation which is independent of n 

zf ~ {e ± pf' z , (31) 



and obtain 



(32) 



In other words, this result means that the scaling relation depends on the strength of the external magnetic field 
with an orientation close to the (xMy) plane for strong Bo, but also on the scales itself with a direction getting 
closer to the (xMy) plane at small scales (small r). This dimensional analysis will be used below to derive the 
unique expression of the vectorial law for anisotropic MHD turbulence since relation (32) gives the following 
dimensional small-scale constraint 

(e ± p) 1/3 ( P \ m 



which leads to n = 4/3. Note that for other types of fluids the value of n may be different [9]. 



5. Exact vectorial law 



Following the critical balance idea we shall rewrite expression (26) for n = 4/3 which gives 

- S (e)e ± re T =F^(r), (34) 
with g{6) = 2cosei(X)/X, 

V3 J \/l + (4/3) 2 tan 2 6> 



and 



I(X) = J y/l + J!/3 dX (36) 





35 7 V J 35 V J 35 V J 



It is the final form of the exact law. We see that the vectorial law has a form close to the isotropic case (4) with 
a scaling linear in r. However, we observe a #-angle dependence which reduces the degree of universality of the 
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law. From an observational point of view this prediction turns out to be interesting since in the solar wind the 
measurements are naturally made at a given angle. Numerical estimate of the function g(9) gives a slight variation 
from 2 to 16/7 for respectively 6 — to 7r/2. It is important to remark that this law is valid for any r and 6 which 
means that we may describe the entire correlation space. Note that the law derived here implies only the mean 
Elsasser energy dissipation rates per unit mass s which makes a difference with other types of universal results 
like in wave turbulence where the spectra may be expressed in terms of directional energy fluxes (like P* or P*) 
[25,26]. 



6. The two-dimensional limit 



It is interesting to analyze the small 6 limit for which the energy flux vector is mainly transverse. For this limit, 
we obtain after a Taylor expansion 

I{X) ~ X + |x 4/3 , (37) 

o 

and then after substitution 

-2 ^1 + |tan 2 6^ e ± pe T ~ F^r) . (38) 

This relation tends asymptotically to the scaling prediction for 2D MHD turbulence which may be obtained 
directly after integration (and application of the Green's flux theorem) of expression (3) over a disk with only 
transverse fluctuations. This result shows in particular how close we are from a two-dimensional turbulence. 



7. Discussion and conclusion 



The interplanetary medium is probably the best example of application of the new exact law in terms of Elsasser 
fields. Indeed, it is a medium permeated by the solar wind, a highly turbulent and anisotropic flow which carries 
the solar magnetic field [27,28] . Several recent works have been devoted to the analysis of low frequencies solar 
wind turbulence in terms of structure functions by using the exact isotropic law [29]. A direct evidence for the 
presence of an inertial energy cascade in the solar wind is claimed but the comparison between data and theory is 
moderately convincing because of the narrowness of the inertial range measured. Some recent improvements have 
been obtained by using a model of the isotropic law where compressible effects are included [30] . Even if the result 
seems to be better the hypothesis of isotropy is a serious default. Other applications of the MHD laws (exact or 
modeled) are also found in order for example to evaluate the local solar wind heating [31] along or transverse to 
the mean magnetic field. 

Direct numerical simulations are very important to check for example the applicability of the universal laws 
discussed in the present paper since there are exact as long as the hypotheses are satisfied. For example, in the 
isotropic case it is interesting to note that the constant has never been checked - only the power law. Therefore, 
we are not yet at the same degree of achievement reached for the four-fifth's law for which the constant has been 
recovered experimentally [32]. Then, for the exact vectorial law derived in this paper it is fundamental to check 
not only the power law dependence (actually, a first analysis at moderate numerical resolution of 256 3 shows a 
relatively good agreement with the scaling prediction) but also - and more importantly - the coefficient g{9) 
which is around 2. Only massive numerical simulations like in [33] will allow to take up this challenge. 

The interplanetary medium is an excellent laboratory to test new ideas in turbulence. In that respect, it would 
be interesting to extend the present work to other invariants like the cross-correlation. Recent works have been 
devoted to this problem where the idea of a dynamic alignment between the velocity and the magnetic field 
fluctuations has emerged [14] but the confrontation with solar wind data is still not totally convincing [34] . Since 
most of astrophysical space plasmas evolve in a medium where a magnetic field is present on the largest scale of 
the system the present law has potentially a lot of other applications. 
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Appendix A. Green's flux theorem 

The appendix is devoted to the Green's flux theorem which may be seen as the two-dimensional version of the 
well-known divergence theorem, ft is also called the Normal form of Green's theorem. Let us consider an oriented 
plane curve C and a plane vector field F defined along C. Then the flux of F across C is the line integral 



/ 



F n< 



(A.l) 



where n is the unit vector normal to the curve C pointing 90 degrees clockwise from the tangent direction of C 
(see Fig. A.f ; left) and d£ is an elementary length of curve C. 





Figure A.l. Left: Oriented plane curve C across which the flux of F is computed. The normal direction is oriented 90 degrees 
clockwise from the tangent direction. Right: Oriented plane curve C that encloses a region S. 

If now C is a curve that encloses a region S counterclockwise (see Fig. A.f; right) and if F is defined in the 
plane (on C and also in S), then we have the relation 



F-nd£ = y J V-FdS\ 

c s 



(A.2) 



which means that the flux of F across a closed integral line is equal to the sum of the divergence of F on the 
surface S. It is the Green's flux theorem. 

A short proof of the Green's flux theorem comes as follows. Let us consider the particular case of a rectangular 
closed curve ABCDA whose orientation defines the x and y directions. On the one hand, one has 



F nd£ = 



c 



c 



dy 
—dx 



a u u a 

= - j F y (x, yi )dx + J F x (x 2 ,y)dy - J F y (x,y 2 )dx + J F x (x 1 ,y)dy 



(F y (x, yi ) - F y (x,y 2 ))dx + J (F x [x 2 , y) - F x {xi, y))dy . 
On the other hand, one has 

%2 V2 

V • F <IS - 

V2 x 2 

[ r i.r-2 .11)- /•:, [ .rj . //) )(/// + J(F y (x, y 2 ) - F y (x, yi))dx , 
yi xi 
which is equal to the flux. 



x 2 V2 

J J V-FdS = j j (d x F x + d v F y )dxdy 



(A.3) 

(A.4) 
(A.5) 
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